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We write the spherical curl transformation for Trkalian fields using differential forms. Then we 
consider Radon transform of these fields. The Radon transform of a Trkalian field satisfies a corre- 
sponding eigenvalue equation on a sphere in transform space. The field can be reconstructed using 
knowledge of the Radon transform on a canonical hemisphere. We consider relation of the Radon 
transformation with Biot-Savart integral operator and discuss its transform introducing Radon-Biot- 
Savart operator. The Radon transform of a Trkalian field is an eigenvector of this operator. We 
also present an Ampere law type relation for these fields. We apply these to Lundquist solution. 
We present a Chandrasekhar-Kendall type solution of the corresponding equation in the transform 
space. Lastly, we focus on the Euclidean topoiogicaiiy massive Abelian gauge theory. The Radon 
transform of an anti-self-dual field is related by antipodal map on this sphere to the transform of the 
' self-dual field obtained by inverting space coordinates. The Lundquist solution provides an example 

04 ' of quantization of topological mass in this context. 

I. INTRODUCTION 

<^ . 

The eigenvectors of the curl operator are called Beltrami fields. These are specifically called Trkalian if the eigenvalue 
] is constant. These arises in different areas ranging from fiuid dynamics and plasma physics to field theories. 
Qh- The field theoretic examples of Trkalian vectors are the Euclidean topoiogicaiiy massive^"— Abelian gauge fields and 
fin ' force-free^ magnetic fields. The topoiogicaiiy massive gauge theories^"— are qualitatively different from Yang-Mills 
-j-J : type gauge theories besides their mathematical elegance and consistency. In this context, Trkalian type solutions on 
^ ' 3-sphere S'^, anti-de Sitter space H'^ and other spaces in connection with contact geometry are discussed in Refs. 

. The Curl transformation, as an interesting tool, is introduced in Ref. The (spherical) curl transformation 
is developed for force-free magnetic fields using vectorsJ^'^*^ The spherical curl transform and Radon transform are 
I [ equivalent descriptions.^*' More precisely, the spherical curl transformation is a Radon probe transformation.^^ Then 
^ ' the spherical curl transformation is applied to the Euclidean topoiogicaiiy massive Abelian gauge theory on MP in 
t~~~ [ Ref. 0- A complex spinor formalism is also developed in Refs. [l2lil3l . 

OO I The Radon transformation and other variants of it have provided valuable insight into problems in different areas 
ranging from tomography to twistor theory. The Radon transformation can be defined for differential forms of certain 
orders in various dimensions j ^ 
' In the next section, we shall first write the spherical curl transform for Trkalian fields using differential forms. Then 
, we shall discuss Radon transform of these fields in its own right. The Radon transform of a Trkalian field satisfies a 
corresponding eigenvalue equation on a sphere in transform space whose radius is determined by the eigenvalue. This 
correspondence can also be inverted. The field can be reconstructed using knowledge of the Radon transform on a 
canonical hemispherei^^ 

In Section Hm we shall consider the relation of the Radon transformation with Biot-Savart^^ (BS) integral operator. 
We shall also present an Ampere law type relation^"^ for the Trkalian fields. Then we shall discuss Radon transform 
' of the BS integral introducing Radon-Biot-Savart (RBS) operator in the transform space. The Radon transform of 
a Trkalian field is an eigenvector of the RBS operator. 

In Section llVl we shall first present applications of our constructions on Lundquist -'^^ (L) solution. Then we shall 
briefly present Chandrasekhar-Kendall^^ (CK) method with circular and elliptic cylindrical solutions. We shall make 
use of an analogous method for finding solutions of the corresponding equation in the transform space. This yields 
an eigenvector of the RBS operator. 

The Trkalian fields yield interesting solutions of not only the topoiogicaiiy massive gauge theories^"^ but gravity^S 
as well. We shall focus on the Euclidean topoiogicaiiy massive Abelian gauge theory in Section [V] The Radon 
transform of an anti-self-dual potential (or field) is related by antipodal map on the sphere to the transform of the 
self-dual potential obtained by inverting space coordinates. 

To the knowledge of the author, L and CK solutions have been overlooked in the context of topoiogicaiiy massive 
gauge theories. The L solution provides an example of quantization of the topological massi"— in this context. 
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We provide the necessary definitions and proofs in an appendix. 



II. THE RADON TRANSFORMATION OF TRKALIAN FIELDS 

A. The spherical curl transformation 

We shall write the spherical curl transformation using differential forms in this subsection. The formulation in 
differential forms follows the same line of reasoning with formally analogous manipulations as given for vectors i^— 

1. The Moses coframe and the eigenbasis 

We use the usual correspondence between vectors and 1-forms: Q\{k) -H> uj°-{k) on R^. We shall also make use 

of the correspondences Vx O *d, V- ^ d*, dot product: (•) O (A*), cross product: (x) ^ *(A) between basic 
operations on vectors and differential forms. 

The Moses^ eigenbasis 1-forms 



Ah-: 



(27r)3/2 



of operator *d 



*dx''{x\k) - Xkx^{x\k), 
d*x''(a;|fc) = 0, a = 1,2, 



d*x^(x\k) 



1 



(27r)3/2 



z/ce*'"-^ * 1, 



form an orthogonal and complete set 



(1) 



(2) 



[x\x\k'),x''{x\k)) = J x''{x\k)A*x\x\k') = 5'^''S{k-k'), 
J2 J \x^ix\k)><x^x'\k)\d^k = I6ix~x'). 



(3) 



Here we use inner product and interior product (bra-ket) notation for 1-forms. 
The complex-valued Moses coframe 



fci(fci + i\k2) 
kik + fca) 



dx^ 



k2{ki + iXk2) 
k(k + ks) 



iX 



dx^ + ^-^±^dx^ 



a = l:A = l,a = 2:A = -l, 



(4) 



LJ^{k) 



-K, k — kidx^ + k2dx^ + k^dx^ , k = |fc|. 



is dual to the basis {Q\{k)} in Fourier space.* The helicity states are given asa = l:A = l,a = 2:A = — 1, 
a = 3 : A = with our conventions. The basis 1-forms of the Moses coframe respectively satisfy the orthogonality 
and completeness relations 



A *uj\k) =< Lj\k) , io" {k) > *1, 
Y,\^^ik)><uj'^ik)\=I, I ={6,,). 



(5) 



The coframe endows with the standard metric 
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ds^ = rjabio^oj'' (6) 




where 



{r^ab) =1-1 I , (g.,) = {5,,). (7) 

We can easily write the foUowing relations 

a;i(fc) = -Zu2(fc), = -ZU^(fc), a;^(fe) = ^^(fc), 

a;"(fc) K = k/k, (8) 

^ _^!l±l^ ^a(^)^ a = 1 : A = 1, a = 2 : A = -1, 

Ki — i\K2 

uj^{k) A uj^ik) = iA * w°(fc), a=l:A=l,a=:2:A = -l, a = 3;A = 0, 

a;°(fe) A Zu''(fc) = iAJ"'' * w^(fc), a, = 1 : A = 1, a, 6 = 2 : A = -1, (9) 

a;^(fe) A *a;''(fe) = 0, a = 1,2, 

u}^{k) A *w'^(fc) = *1, 

a;^(K) A A = * 1. (10) 

Here (~) denotes complex conjugation, (*) is the Hodge dual, *1 = dx^ Adx^ f\dx^ , uj = (/«)| = i, \/\rf\ — i, = 1 
and ** = 1. We refer the reader to Refs. i8i- tl0i.2lM22. for a comparison with formulation in vectors. 

2. Trkalian field 

The curl transform of a 1-form field is given by an expansion in terms of these eigenforms as in the vectorial case. 
We can think of this as a plane wave expansion in the coframe 
A Trkalian field *F 



can be expressed as 



:d*F-z/*F = 0, (11) 



^Fix)= *"^(^): (12) 



where 



* ''Fix) = - [ x''{x\k)fa{k)d^k, (13) 
9 J 

excluding the divergenceful component i^i^iiiS The factor 1/g is introduced for the sake of a proper strength for the 
gauge potential in topologically massive gauge theory.^ This can be taken as 1 for general Trkalian fields. Then the 
curl transform of *°F(cc) is given as 



(14) 
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We find 



faik) = ^^^p^-.(fc), (15) 

using a radial delta function.— An arbitrary solution is given entirely in terms of its transform on the sphere of radius 
k = Xv = \v\. Furthermore, only the eigenforms for which A = sgn{v) contribute to the field (jl2p .^ The expansion 
simplifies into 

*°-F(x) = - / x''{x\>^vi^)sa{XvK)dn (16) 
9 J 

1 1 



(27r)3/2 g 



e"""'''L0''{K)Sa{XvK)dVL, 



where d^l is the spherical area element and k = fc/fc is a unit vector in transform space. Therefore a solution can be 
defined entirely by the value of its curl transform on the unit sphere in transform spacei^ We call Sa the spherical curl 
transforn*i^ in order to distinguish it from the full curl transform fa ■ 
The spherical curl transformiS is given as 

where 

F^{p,k)^ y"*"F(a?) A*w"(K)(5(p-K-a;) (18) 

= < w°(k), J *°F(a;) > S{p- K-x)*l, 

is the Radon transform of the component Fa{x) =< *'^F{x) > of *°'F{x). The Radon transform of a function 

(see Appendix lA ip is basically defined as the integral of that function over the plane at a distance p = k • a; to the 
origin with unit normal k. 
We find 

F;^{p,k) = (27r)i/2 -\< oj^iK), \e'^"Puj''{K)sa{XiyK) + e-'^''Puj''{-K)sa{-XiyK)] >, 
g ^ ■' 

substituting (fT6)) in (fT8|) and following the reasoning in Ref [l3 with our notation. We shall not repeat these 
calculations here since they are formally analogous. Then we construct 

* ''F'^ip, k) = (2^)1/2 - \ [e'^''Pu:''{K)sa{XvK) + e-'^''Puj''{^K)sa{-XvK)\ , (19) 
9 

using the completeness relation ([S])- We can write the spherical curl transform p7)) as 

^a{X>^'^) - jArnd'^^^'"^'''' * [*''F'^{P, «) A *i0'^{K)] (20) 

(iTTj-^/^ 

- (^^•9^'^"^"'^' < ^"('^),*"^''(P,'^) > ■ 

This is a Radon probe transformatiooi^ (see Appendix lA ip . The Radon transform *°'F^ {p, k) is composed of both 
helicity components of the field This satisfies the relations 

^^F'^{-P,-k) = *''F'^{p,k), *''F'^{-p,k) ^*''F'^{p,-k). (21) 
See Appendix [B] for equivalence of the spherical curl and Radon transformations. 

We have given a formulation of the spherical curl transformation in terms of differential forms. We shall use 
vector formulation below which is more suitable for Radon transformation in three dimensions. One can use the 
correspondence between differential forms and vectors for translation. 
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B. The Radon transformation 

The Radon transform F^{k • a;, /«), (|T9| ioi p — k - x which is the integral of the Trkahan field over the hyperplane 
orthogonal to k containing x satisfies 

V X F'^{K-x,K)-iyFj^{K-x,K) =0, (22) 

([TT|) . where V • F^{k ■x,k) = 0. 

We shall prove that the Radon transform of a Trkalian field satisfies a corresponding eigenvalue equation in the 
transform space. We shall interchangeably use [ ]^ or ??.[ ] for denoting the Radon transform. (See Appendix I A 11 ) 

Definition 1: The operator T is defined as 

d _ d 
op op 

Here d/dp is called the (infinitesimal) parallel displacement operator of the planej^ 
Proposition 1: The Radon transform intertwines: 

a) the curl operator Vx and Fx 

7^[v X v{x)]{p, k) = r X 'ri[v{x)]{p, k), 

h) the divergence operator V- and T- 

n[v ■ v{x)]ip, k) = r • n[v{x)]ip, «), 

c ) the gradient operator V and T 



Proof: We have 



n[Vf{x)]{p,K)^Tn[f{x)]{p,K) 



^^5{p~K-x) = -K,5'{p-K-x), ^ = 



s = p — K ■ x, (23) 



d Q Q s 

6'{p- K ■ x) = —6{s) ^ —S{p- K- x), ^ = 1' 



which yield 



VS{p- K ■ x) = -kS'{p - K- x) (24) 

— —K-^S(p — K ■ x). 

op 

Thus2i^ we find 

n^J{^m z^) = ^5(P - ■ x)d'x (25) 

= K,^n[fix)]{p, k). 

The domain D of integration is the whole space: D = M.^ unless otherwise stated. This leads to componentwise proofs 
of a, b and c. 
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a) More precisely, consider a S-dimesional (finite) region D bounded by a surface S = dD. We have 

/ V xV{x)S{p- K-x)d^x^ I V x[V{x)S{p-K-x)]d^x- VS{p-K-x)xV{x)d^x (26) 
Jd Jd Jd 

= - / S{p - K ■ x)V{x) X ds - VS{p - K ■ x) X V{x)d^x. 
Js Jd 

The first integral vanishes if we consider functions V{x) which vanish at infinity, as the region D is extended to whole 
space. (For a finite region we could use vectors normal to the surface.) Thus 

V X V{x)6{p -K-x)^ - V(5(p -K-x)x V{x), (27) 

(under integral). Hence we find 

7l[V xV{x)]{p,k) = I V xV{x)S{p- K-x)d^x (28) 
Jd 

^Txn[V{x)]{p,n), 



using ^ and ((24)) . 

The proofs of b and c are based on a similar reasoning. The ( X ) product is respectively replaced by (•) and ordinary 

Saduct, with similar boundary conditions at infinity (for a finite region, with certain boundary conditions). See Ref. 
for alternative proofs with different techniques. □ 
Definition 2: The operator is defined as ~T ■ T. 
Proposition 2: The operator acts on f^{p,K) as 



T^f'^{p,K)^{dydp')f'^{p,K). 



Proof: We simply have 



rV«(p,K) = r-|r/«(p,K)l 

= ^/^(P,'«), K-K^l. 

□ 

This operator corresponds to □ = which is intertwined^^ with Laplacian V^. This acts componentwise on vectors 
V^{p, k) analogous to the action of on V{x). 

Proposition 3: The operator T satisfies the following identities: 

a) Tx [r/^(p,K)] =0, 

b) r- [r X v'^{p,k)]^o, 

cj r X [r X v^(p, k)] ^ T[T ■ v^{p, k)] - fV^(p, «). 

Proof: These follow as 

a) r X [r/^(p, H,)] = X ['^^/'^b, «)] 

-d-p^'^'^'^w (^''^^ 

= 0, K X K = 0, 
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b) r • [r X k)] = ■ [k-^ X y^(p, k)] 

= 0, X V^(p,k)] = 0, 

c) r X [r X «)] = x [k^ x «)] 

= r[r.y^(p,K)]-r2y^(p,/^). 

□ 

Proposition 4: T/ie Radon transform intertwines the eigenvalue operators (Vx) — v = and (Fx) — = for 
constant eigenvalues. 

Proof: The proof follows from the linearity of the Radon transformation. □ 
We can easily check the Radon transform (fT9|) satisfies 



^{[k-V'^{p,k)]k-{k-k)V'^(p,k)} 



r X [p, k) - ,yFl^ [p, (29) 

[in differential forms: dp = k = —uj^, ^ sincep — k-x]. We also have: K-F^{p, k) = which leads to r-F^{p, k) — 0. 
We can write (f29l) as 

(p, K) + ,yKX F^{p, k) = 0. (30) 

A change of sign: p — > —p or k — > —k (but not both) leads to a change of sign in the eigen- value: v — > — in 
(l29t . We shall discuss this again in Section |Vl The Radon transform T U^{p, k) of an arbitrary gauge transformation 
VU{x) is normal*^ to S^. 

We need the adjoint Radon transformation for discussion in the reverse direction. 

Definition 3: The adjoint^^ Radon transform of a function G{p, k) on the transform space is 



G^\x) ='R)[G{p,k)]{x) = [ G{K-x,K)dn. 



The function G{k ■ x, k) is to be taken as G{k ■ x,k) ~ F''^{k ■ x, k). The Radon transform TZ integrates over the 
set of points in a hyperplane, the adjoint transform 7?J integrates over the set of hyperplanes through a point.— 
Proposition 5: The adjoint Radon transform intertwines: 

a) the operators Tx and curl Vx 

7e.t[r X G{p,k)]{x) = V X 7^t[G(p,K)](a;), 

b) the operators T- and divergence V- 

n^[r-G{p,K)]{x) = v- -r) [g(p, k)]{x), 

c) the operators T and gradient V 



lV[T<ty{p,K)]{x) = VTV[(^{p,K)]{x). 
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Proof: The componentwise proofs of a, b and c follow from the identity 



_d_ 



-R) k)] {x) = ^ J^J{K-x,K)dn (31) 

— [ K^-^(h(p, K)dfl, p—K-X 



= n^^,^-^<Pip,K)]{x). 



□ 



The adjoint Radon transform also intertwines^S- the operators and V^. 

Proposition 6: The adjoint Radon transform intertwines the eigenvalue operators (Fx) — = and (Vx) — = 
for constant eigenvalues. 

Proof: The proof follows from the linearity of the adjoint Radon transformation. □ 
Thus, we have a correspondence between the (constant) eigenvalue equations (V X ) — = in the physical space 
and (Fx) — = in the transform space. Note that we should employ (— l/87r^)V^ on 7iJ[F''^{p, k)] for a complete 
inversion. This simplifies for a Trkalian field. Further one should consider the appropriate class of functions for the 
Radon transform and its adjointi^ (See Appendix lA 11 ) 



1. A refinement 

The Radon transform (fT9l) of a Trkalian field is composed of both helicity components of the field. The inverse 
Radon transform (jA4p requires knowledge of the transform evaluated over the entire sphere,^- However the knowledge 
of a helicity component evaluated over the entire sphere suffices to reproduce the knowledge of the other component. 
Thus construction of the field can be achieved using the knowledge of a single component evaluated over the entire 
sphere as in ^TE\i . Equivalently, the field can be constructed using knowledge of both components, i. e. the Radon 
transform (|19p . on a canonical hemisphere which can consist of disconnected parts. For scalar fields, a similar result 
is proved in Ref. [l^. Here, this easily follows as one considers both helicity components in the Radon transform of a 
Trkalian field. This also removes redundancy of double-covering of spacei^^ 

A canonical hemisphere if is a (Lebesgue measurable) subset of the unit 2-sphere whose area is 27r such that 
if the tip of a unit vector k based at the origin is in the hemisphere then the tip of — k is not in the hemispereJ^ 
The complementary set H' in is a canonical hemisphere too. A canonical hemisphere can consists of disconnected 
parts. 

If we substitute with p — k ■ x in (|A4p , we find 



^^(^) = 'tAtH - I [e*^"" "QA(K)sA(Az.K) + e-*^'"^-^QA(-K)sA(-Ai./«)] d^ (32) 

(27r)-^/^ g Jh 

= lAzH - I [e*'"" "QA(A^)sa(Ai^/^) + e-'''-'^ '-Q\{-K)sx{-\i^K)\ dQ, 
(27r)'^/^ g Jh 



as K — > K. in H', using the identities ([2T|) . This can also be inferred from ([161) . 

We can see this on a simple example choosing sa(Az^k) = {2TT)^^'^gS{K. — Kq). We decompose the sphere into 
two complementary canonical hemispheres H and H' where the tip of Kq is contained in H: Kq e H. If we take 
sa(— Ai^k) — (27r)'^/^g(5(/« — k'q), (k — > —k: Hq — > k'q — —Kq) then k'q e H'. We find 



Fx{x) = e'^'^'^-^QxiK^) (33) 

using (|32p or equivalently inverting it and integrating over H' . 

It is possible to define a left and right inverse of this refined (inverse) transformation if it is applied to Radon 
transforms of functions (see Appendix lA 2|) . Alternatively, one can introduce inverse of this refined transformation 
redefining the Radon transform and the set {F^(p, n)} to require that k be restricted into a canonical hemisphere. 
Then 7^^^ is both a left and right inverse of 7^J^: 7^^'^7^J^[F](a;) = F{x), IZhT^'hIF'^Mp, «) = F'^{p, i^)'^ 
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III. RADON TRANSFORM AND BIOT-SAVART INTEGRAL 

In this section, we shall consider the relation of the Radon transformation with Riesz potential and BS integral. 
We shall also present an Ampere law type relation for Trkalian fields. Then we shall discuss Radon transform of the 
BS integral. 

We can easily^^""— show 

■R.^n[F]{x) = 8Tr^I^[F]{x). (34) 
The integral on the right-hand side which is defined as the convolution of F with Riesz kernel K2 — (l/47r)|a;|^^ 



I^[F]{x)^iK2*F){x) (35) 



1 f F{v) ,3. 



47r Jjj \x ~ y\ 

is called the Riesz potential for the vector field Fix). We remind: D — for the Radon transformation. The 
operator — is an (left) inverse of the Riesz integrali^ 
We find 

X ■R)n\F]{x) = ^■K^BS[F]{x), (36) 

where 

BS[F]{x)=V.,xI^[F]{x) (37) 
= f / F{y) X ^-^d^y, 

is the BS integral operator.-ii The vector potential for BS[F]{x) is given by the Riesz potential (|35|) .— 
We immediately see that 

• X 'R)n[F]{x) = Svr^V^ • BS[F]{x) = 0. (38) 

We also find 

X X n'^n[F]{x) = Stt^V^ X BS[F]{x) = 8tt'^F{x), (39) 

if F is divergence-free and it also vanishes at infinity (for a finite region D, it is tangent to the surface S — dD 
bounding this region). (See Appendix [Cl) Hence the BS integral is divergence-free and the curl operator is a left 
inverse of this under these conditions 

We can associate an eigenvalue equation for BS operator (the eigenvalue being reciprocal of i^) with the Trkalian 
fields only in this caseJ^!^ In general, the curl can be taken to be the inverse of an appropriately modified BS 
operator J^S^, 



As for the Trkalian fields, we find 



■R^1Z[Fx]ix) ^ 8n^^Fx{x) 



(40) 



using Definition [3] and (IT9|) . (fT6|) . This can be inferred from (jA4p using which is intertwined with V'^. Note that 
the conditions in (p9| are satisfied for a Trkalian field if its Radon transform exists (see Appendix lA ip . The equation 
(jiO]) yields 

V, X V, X TZ^TZ[Fx]{x) = Sn^Fxix). (41) 



K. Saygili: Trkalian Gelds and Radon transformation 



10 



This reduces to the inverse Radon transform (jA4l) using (l22l) on the left-hand side. 
We find an Ampere law type relation^i^ 

$ = lyQ, (42) 
integrating the equation (ITT|) over a surface S bounded by the curve C — dS. Here 

Q = [ F[x) -ds, $ = / [V X F{x)] ■ ds, (43) 

Js Js 

are respectively the flux of F{x) and V X F{x) through the surface S. The flux $ reduces to 



$ = / F{x) ■ dl, (44) 
Jc 



the circulation of F{x) on C. Thus the flux of a Trkalian field through a surface is determined by its circulation on 
the boundary of the surfacej^i^ Note that the equation (|22|) yields an analogous relation for F^[k ■ x,k). 

A. Radon transform of the Biot-Savart integral 

We shall make use of Fourier slice theore m^^'^^ 

T[F'^{p, K)]{k, k) = 2TTT[F{x)]{kK), (45) 

for finding Radon transform of the BS integral. On the left-hand side stands for 1-dimensional Fourier transform 
in the first variable whereas on the right-hand side it is a vector Fourier transform in 3 dimensions. We also need 
Fourier transform of the Ries a^^i^^ integral 



1 

k 

1 1 
2^k 



ni'[F]ixmk) = ^JP[F(a;)](fc/.) (46) 
\T[F'^{q,K)]{k,K). 



Here the last line follows from (|45 
We find 



:F{/2[F]^(p,^)}(fc,«) = 2nJ-{l'[F]ix)}ikK,) (47) 
= ^JP[F^(g,K)](fc,,.), 



using and This yields Radon transform of the Riesz integral 



'R,{I^[F]{x)}{p,K) = T-'{^T[F'^{q,K)]{k,K)}{p,K) (48) 



1 

2^ 



F"iq,K)dq. 



This equivalently corresponds to 7l-7?JlZ[F]{x), ([34]).— This is a convolution for the Radon transformation<2^i21 The 

operator — F^ is an (left) inverse of the Radon transform of the Riesz integral. We define Radon-Biot-Savart integral 
operator as 

RBS[F-^{q, K)]{p, k) = BS[F]'^{p, k) (49) 
= X J^-\^T[F'^{q,K)]{k,K)}{p,K) 

= J^-'{^J^[T, X F'^{q,K)]{k,K)}{p,K), 



K. Saygili: Trkalian fields and Radon transformation 11 



the Radon transform of the BS integral using Proposition [T] and (|37)) . (l48l) . 
We immediately see that 



r-RBS[F'^{q,K)]{p,H!,)=0, (50) 

using Proposition |3l We find 

r X RBS[F^iq, k)]{p, k) = -k, X [k X F'^{p, k,)] (51) 

= -{[k ■ F^(p, k)]k - [k ■ K)F'^ip, k)}. 

This yields 



r X RBSiF'^iq, k)]{p, k) = F^(p, k), (52) 

if K ■ = 0. This condition leads to F • — [as can be inferred from (j49|) ] which corresponds to V • = in 
the physical space. Thus we conclude that Fx is a left inverse of the RBS operator, if k ■ f ^ = 0. We can associate 
an eigenvalue equation for RBS operator with the Radon transform of Trkalian fields. 
The RBS operator also acts as a left inverse 



RBS[r, X F'^{q,K)]{p,K) = F'^{p,K), (53) 

of the operator Fx for those vector fields F'^ with F • F'^ — 0. For example, the Radon transform ([T5|) of a Trkalian 
field is an eigenvector 



RBS[F^{q, K)]{p, k) - -Ff (p, k), (54) 

of the RBS operator with reciprocal eigenvalue. The kernel of the RBS operator consists of those vector fields for 
which F X F^ = i. e. gauge terms FJ7^, (Proposition corresponding to gradient fields VC/ in the physical 
space 



IV. EXAMPLES 



In this section we shall first present applications of our constructions on the L solution. Then we shall briefly 
discuss the CK method with circular and elliptic cylindrical solutions. We shall make use of an analogous method for 
finding solution of the corresponding equation in the transform space. 



A. Lundquist solution 

The Li^ solution is 



Fl = Fo[Ji{iyr)eg + Jo{iyr)e,], (55) 

where J™ is the Bessel function of order m and Fq is the strength of the field on z-axis. Note this is of helicity X = iJ^ 
The spherical curl transform p7|) of ([55)1 is given by 



si{vk) = -2^/2(2^)1/2^^^ S{K,)e-'^. (56) 

The radial degree of freedom k is also taken into accounti^S Here: k = kr cosT/^ej; + kr smipBy + k^e^ in the circular 
cylindrical coordinates in transform space. Thus the spherical curl transform of the L field lives on the equatorial 
circle of the sphere,—- 
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This solution provides a simple example of Ampere law. Consider a disc S which is bounded by the circle C = dS 
of radius R, centered at the origin in xy-j>\ane. We find 



$ = i^Q = 2ttFo RJi {vR) , (57) 



32]). Note that $ = Q = at zeros of Ji{vR). 
The L field is an eigenvector of the BS operator 



BS[Fl\{x) = ^FUx), (58) 



37)) . This calculation is given in Appendix iDl 
We find 



F^(p, n) - 2mFo^S{t,,)ie"'PL + e-'^^L'), (59) 

where L — sin^Sx — cosipBy — ie^ and L' ^ — sin ^/'e^; + cos^By — ie^- Note that we have taken the cases: k negative 
or positive into account here. We immediately see that F]^{p, k) satisfies the equation ([29]) anticipating the effect of 
6{kz)-^'^ This is also an eigenvector (|M)) of the RBS operator. 



B. Chandrasekhar-Kendall cylindrical solutions 



Chandrasekhar-Kendalli^ developed a method for deriving solutions to V X = vF from the scalar Helmholtz 
equation. See also Refs. Issl - lssl 
We find 



(V2+i.2)F = 0, (60) 

applying curl on'V xF — vF = 0. Thus every solution ofVxi^ — i^i^ = Oisa solution of (pO)) . although the 
converse is not necessarily true. If ^ is a scalar function satisfying the Helmholtz equation 

(V2 + i^2)* = 0, (61) 
then three independent solutions'^- of the equation ([60l) are 



p = Vx(*a;), q = -V X p = -V X V X (*u;), r = V*. (62) 

Here w is a fixed vector (of unit norm) . The vectors of type p and q are respectively called toroidal and poloidal.— 
We also have: V X <7 = vp. Hence we find: V X (p + </) = v(j) + q). Therefore, the most generaF-^^ solution of 
V y. F — V F — Q among the solutions of (160]) is 



F = p + q = V X (*a;) + -V X V X (fa;). (63) 

The scalar ^ is called Debye potential. 3- This method has been frequently used for discussing circular cylindrical, 
spherical and other type solutions.^ For the relation with Wu-Yang type solutions in the topologically massive Abelian 
gauge theory see Refs. 00. 

The CKi^ solution in cylindrical coordinates is given by 

F^-[(tV X (*e;,) + V X V X (^e^)], (64) 
where V X F — ctF = 0. In the circula r^°i^^ cylindrical coordinates 
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^{r, e, z) ^ Jr,^{vr)e""^-'^\ a'^ = v'^ + e. (65) 

This reduces to the L solution ([55]) for to = 0, /c = 0, {Fq — —v^). In the elhptic cyhndrical coordinates the solution 
is given by (1641) with 

^{u,v,z) = U{u)V{v)e-^^'-, = k'^ = v^. (66) 

Here U{u) and V{v) are respectively the modified Mathieu function and the Mathieu function. For k = Q^a = m = v 
and this solution reduces to Vandas-Romashets^^il^ (VR) solution which is an elliptic generalization of L. 

1. Chandrasekhar- Kendall method in transform space 

We can adapt the CK method in the transform space for finding solutions to F X G = G. Briefly, the solution of 

T xG{p,k)-vG{p,k) = Q, (67) 

that is to be found among the solutions of 

{V"^ + v^)G{p,k) = Q, (68) 

is 

G(p, k) = r X (fa;) + ir X r X {^ijj). (69) 

Here ^'(p, k) satisfy 

(r^ + i/2)5r ^ 0, (70) 

and a; is a fixed vector in the transform space. We have ignored the superscript R for Radon transform in '^uj. We 
remark a few important issues. First, the Radon transform is defined on X R [or {M? — {0}) X M], (see appendix 
IXTj) . Secondly, G{p,k) may consist of distributions. The third issue is imposing boundary conditions in the transform 
space. 

The first term in G{p,k) = T X (fa;) with * = [{2tt)^ /i^^][e'^''PS{K- Ko) + e-^^''P6{K + Kq)] and uj = Qa(ko) 
immediately yields 



G{p, k) = (27r)2 _ [e'^-PSiK - Kq) + e''^''P5{K + Ko}] Qa(a^o), (71) 

which satisfies (|67l) . This is the Radon transform: G{p,n) = F'^{p,K,) of the solution F{x) = e^^'"^°''^Qx{Ko).^'i^ 

As another example, consider the Radon transform ([59|) of the L field. We can reproduce G(p, k) — F]^{p, n) from 
the first term in G{p, /«) = T X (^-i^i + ^2(^2) choosing *i = {2TriFn/v^)S{K^)e'''P , = {2niFo/iy^)S(K:,)e-"'P 
and u)i = L, u)2 — L'. The Radon transform of Jo{i'r)ez leads to a simpler choice: 'I' — {2TTFo/i/^)S{Kz){e^'^P + 6^^'^^), 

UJ = Bz. 

We can easily show that G{p, k), ([M)) is an eigenvector 



RBS[G{q, k)]{p, k) = -kx[hx H{p, k)] = ^G{p, k), (72) 

of the RBS operator (|49)) . if a; is a fixed vector: dpU} = 0. Here H{p,k) — 'i'uf + {l/i')T X (^t'o;) is a potential for 
G{p, k): G{p, k) = r X H{p, k) and we use (f70| . 
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The CK method also enables us to write a simple integral expression for the solution of (|67)) . We employ in ([69)1 
integral expressions for simple harmonic oscillations (1701) . If we assume ^'^(p, (fc = 1, 2) are separable in p and k, 
these are given as 



Uk[K) 



(73) 



on the infinite interval: —oo < p < oo, using Laplace kernel: e^^ in the complex (-p\a.ne.^*'*'^ Here are arbitrary 
loops around the poles C = ±iXi^- This is basically the Cauchy integral representation formula for e'^^^'^P. We find 



G(p, k) 



47ri 



{iXK X u)i — K X K X ivi) / 7—d( + {—iXn X uJ2 ~ k X k X 1^2) I 

Jc- C - Jc 



(74) 



using ((73)) in ([69]) and rearranging it. This is also oscillatory in p. We have included the k dependent factors Uk{K) 
in the arbitrary functions u)k{K,). These are decoupled from the integral expressions. 

For example, if we choose a;i(/«) = [(2Tr)^/'^ /{giy'^)]Qx{K)sx{XiyK) and iV2in) — [{2?:)^^'^ / {gv'^)]Q x{~k)sx{~Xi^k,) 
in accordance with the Moses frame, we find G{p, k) = F^{p, k), (IT51) . 

We can reproduce the L field ([59]) choosing uji = S{Kz)ez and u}2 — S{Kz)ez- We can also use uji = S{Kz)e^ and 

U}2 = 5{Kz)ea, {Kz = cosa, K — Br) On S^. 

If we substitute (1741) in the inversion formula (jA4p . we find 



Ffa;) 



1 



327r3i 



/ - — [ {iXn X LJi - K X K X uJi)eS'^''"d^dC (75) 

+ / - — 5—— / {-iXk X u>2 - X K X u)2)e'^'^''"dnd( , 
Jet C + Js2 



using the operator which is intertwined with V^, (|68p and interchanging the order of integrations. 
If we choose 

= -i{2n f-^ [aS{K, - ki)Ei + b5{K. - K2)E2 + c5{k - kz)Ez\ , 
UJ2 = -i{2TTf\ [a5{K + K,i)Ei + bS{K + K2)E2 + c5(/« + /«3)-E3] , 



(76) 



where a, b, c are arbitrary constants, iXhii X Ei — Ei and ■ Ei — 0, the equation (j75p yields 



F{x) = -i (ae'-^'^'^i-^£;i + be'^"'''' '^ E2 + ce'^"''' '^ E3) . (77) 

The real part of this reduces to abc fielc^ if we choose: ki ~ k2 = — Xv, Ki ^ Bz, k,2 — e^, — By and 
£^1 = ej; + iXey, E2 = By + iXbz, E3 ^ Bz + iXb^- 

We respectively refer the reader to Refs. I46ll47l for Fourier integral expression on §^ of the general solution to the 
scalar Helmholtz equation and to Refs. |48| - |5C)I for integral representations of the Beltrami fields in physical space. 



V. TOPOLOGICALLY MASSIVE ABELIAN GAUGE FIELD 



The Euclidean topologically massive Abelian gauge field^"'^ is a Trkalian field. ^ Another example of Trkalian fields 
is the force-free magnetic fieldji The Trkalian fields yield solutions of both the topologically massive gauge theories^"— 
and gravity)^ The spherical curl and Radon transform analysis above arc also valid for the Euclidean topologically 
massive Abelian gauge fields on R"^. The gauge potential of the Trkalian field is also Trkalian. Meanwhile a gauge 
transformation corresponds to a vector normal to the sphere. The Radon transform of an anti- self- dual potential (or 
field) is related by antipodal map on the sphere to the transform of the self-dual potential obtained by inverting space 
coordinates. 

The Trkalian fields L and CK are also solutions of this theory. Furthermore, the L solution provides a simple 
example for quantization of the topological mass in this context. 
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A. The gauge potential 

We can easily derive the spherical curl and Radon transforms of the potential following the saniei^, reasoning above. 
We find 



A\(x) — - — I XA(a5|Ai/K)sA(Az/K)df2. 
g V 



(78) 



This satisfies the self-duality relation: Fx — vA\ where F^ — '^y. A\, ([TB|. We also find 

Aj(p, k) = (2^)1/2 i -1 {e'^''^Qx{K)sx{\vK) + e-'^''-PQx{,'K)sx{-\vK)\ . (79) 
These only differ by a factor^ of Xjv from those expressions for F\ . Hence the potential Aj(p, k), (|75|) satisfies 



Ff{v,K)-vAl{v,K)=(), 



where FjJ^ = T X Aj, ^ 
A gauge transformation 



(80) 



A' = A - - VC/, 
5 



of the potentials is given by a curl-free vector 



This yields the gauge function 



VV{x) = i y xo(a;|fc)/o(fe)d3fc. 



(81) 



(82) 



^(.)=(^.i/e'-/o(fc)i.3fc. 



(83) 



B. The anti-self-dual case 



The anti- self- dual case 



V X F + i^F = 0, 



(84) 



of equation (jlip contains an extra factor of (— ) in v. This terminology is motivated by the interchange symmetry 
F ^ vA. This leads to a flip of sign in the spherical curl and Radon transforms of the self-dual field (|lll) . In this 
case, an arbitrary solution is given in terms of its transform on the sphere of radius fc = —\v = \v\. Furthermore, 
only the eigenfunctions for which A = —sgn[v) contribute to the field. We find 



Fx (a;) = - [ Xx{x\ - Xi^K)sx{-Xi'K)dn 
9 J 

''Q\{K,)sx{-XiyK,)dQ, 



(85) 



1 1 

(2^ g 



and 



(86) 
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The Radon transform 

(p, k) = (27r)i/2 i 1 [e'^''PQ^{-K)sx{XiyK) + e-'^''PQ>,{K)sx{-Xiyn)] , (87) 

of the field satisfies: T X — —vF^. Noting Xv ~ > for the self-dual case and —Xv = > for the 
anti-self-dual case, the equation ([57)1 coincides with (|19l) . The difference of the two cases is the relative sign of p and 
K. We shall present a comparison of the two cases below. 
The gauge potential 



Ax{x) = - y x\{x\ - XvK)sx(~Xvhi)dQ., 



for the field ((85|) satisfies the anti-self-duality relation: F\ = —vAx. Its Radon transform 

A^{p, /^) - -(27r)i/2 - 1 [e'^''PQx{~K)sx{XvK) + e'''^''PQx{K)sx{-XvK)\ , (89) 
g 

satisfies: F'^ = -i^A'^. 

C. The Radon transform and duality classes 

We can compare the Radon transforms of self-dual and anti-self-dual fields using a correspondence between them. 
A self-dual field becomes anti-self-dual, or vice versa under inversion of space coordinates, if the topological mass 
is held fixedi^ (Hence the topological mass is actually a pseudoscalar.) Therefore, given a self-dual-field F{x): 
'Vx X F{x) — V F{x), we can write the corresponding anti-self-dual field as F'{x) — F{~x): 'Vx X F'{x) = —v F'{x) 
inverting the coordinates: x — > x' = —x. 

If T is (matrix of) a nonsingular linear transformation: T^^x = x',— then 

P«(p, K) = K|F'(i)l(p. = ■KlF(T-^x)]lj,. k) (90) 

= |r| J F{x')S{p - T^K ■ x')d'^x' 
= \T\F^{p,T^K). 

The inversion of coordinates is given by T = —1, |T| — \det\T]\ = 1 which is an orthogonal matrix: TT^ = 1. Thus, 
the Radon transform of an anti-self-dual field 

F'^(p,/^)-F^(p,-/^), (91) 

is related by antipodal map on the sphere to the transform of the self-dual field obtained by inverting the coordinates. 

In the Moses coframe, we can find the Radon transform of an anti-self-dual field similarly inverting the coordinates 
in ([TSl) . This leads to the inversion n — > — /« in ([T9l) . A flip of sign ([21]) : p — >■ —p or k — > —k (but not both) 
brings * "F^lp, -k) back to with Xiy = > or (gT]) with -Xi^ = > 0. Note that the equation ([IT]) together 
with ^T9\\ yields a representation different from ((87)) for anti-self-dual fields. These coincide upon both: ly — > —v and 

K > —K. 

The Radon transform of an anti-sclf-dual field or potential satisfy the equations (P5)) or (|80p with an opposite 
sign for the topological mass as noted above. A simple example is given by F{x) = e'^'^'"'"-='Qx{Ka) where ^ = ±1 
correspond to the self-dual and anti-self-dual cases. We find 

F^(p, K,) = (2^)2 1 [e'^-PS{K - /iKo) + e-'^''PS{K + ^i^,o)] Qa(«o), (92) 
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using Xv > in the self-dual case (ffTj) and —Xiy > in the anti-self-dual case. This yields T X = ^ivF^ 
anticipating the effect of 5{n =F /i/to)-^" 

Noting the inversion of the refined transformation for a self-dual- field F^'{p,k) on a canonical hemisphere in 
Subsection IIIB 11 one can require k to be restricted into the complementary hemisphere for the anti-self-dual field 
F''^{p,k). (Also see Appendix |XT]) 

D. The Lundquist Solution 

The constant eigenvalue v, called the topological mass in this context, is intimately related to boundary conditions 
imposed on the fiel d^^'^^i'^^ and the global topology determined by them. The quantization of the topological mass 
here arises as a result of a well defined gauge transformation as in the previous^i^ examples. 

The eigenvalue in L solution determines the twist per unit length of the field lines on the z-axis, see for example 
Refs. [5M53I In the topologically massive case, the strength of the field on the z-axis is proportional to the square of 
the eigenvalue. Meanwhile the gauge potential is of strength proportional to this twisting. This analogy suggests a 
natural geometric interpretation for the Lundquist field in this context. 

We find the potential 

A = -FL--Foe,, (93) 

V V 

integrating^ the components of i^L, (|55|) . This satisfies: V X A — vA = F^e^ with an extra term. 

In the topologically massive Abelian gauge theory^^i^ we can make this term vanish: F'j^ — vA! = by a gauge 
transformation: A! — A — {i/g)'^\\\U identifying it with a gauge term: — (j/g)Vlnf7 — {l/v)F[)e^ . This yields 
U = e*"^ which takes values in the group C/(l), choosing Fq = /g- The strength of the gauge potential is given by 
the gauge coupling constant v/g = ng, li v ^ ng^. This leads us to adopt a fundamental scaled of length I — ^-n j g^ . 
We can write v = ng^ d& v — 2TTn/l. If [/ = e*^^'^"/')^ is a single-valued function of z with the fundamental scale 

then n has to be an integer. The fundamental length scale I is the least common multiple of intervals over which 
the gauge function is single-valued and periodic for any integer n, in addition to the fact that it has a smaller period 
V = l/n^ The discussion of physical and topological aspects of this solution and the effect of gauge transformations 
on these would be distracting us from our purposes here. 

The anti-self-dual solution is given by /t — > — k in ((59)l . Also note that the L field defines a contact structure since 
Jo and Ji have no common zerosi^ 

VI. CONCLUSION 

The spherical curl transformation for Trkalian fields is a Radon probe transformation in the Moses frame. We have 
written this using differential forms. This is an example of an integral transform which can be formally written this 
way. 

The point of view here suggests an approach for studying the Trkalian fields in the transform space. The transform 
space representation has certain advantages. First of all, the equation in the transform space provides a concrete 
frame for studying these fields, for example using the spherical curl transform which is based on Moses eigenfunctions. 
This further offers a simplification since the differential and algebraic aspects are separable in the transform space. 
This also exhibits the role of the Moses eigenbasis which yields a helicity decomposition of the Radon transform. 
Meanwhile the field can be reconstructed using knowledge of the transform on a canonical hemisphere. Briefiy, this 
equation provides a simple geometric frame in the transform space for Trkalian fields with a new insight into their 
structure besides practical advantages. 

We have discussed the connection of the Radon transformation with the BS integral. We have introduced the 
RBS operator in the transform space. This is simply given in terms of 1-dimensional Fourier transforms. We can also 
associate an eigenvalue equation for RBS operator with the Trkalian fields in transform space. The Radon transform 
of a Trkalian field is an eigenvector of this operator. Its kernel consists of those vector fields which are normal to the 
sphere. We have also presented an Ampere law type relation for these fields. 

We have presented applications of these on the L solution. Then we have discussed the CK method with circular 
and elliptic solutions. The elliptic one reduces to VR solution. The L solution also defines a contact structure. 

We have presented CK solution of the corresponding eigenvalue equation in the transform space. The CK method 
in transform space also leads to a simplification eliminating the differential operations and hence reducing the solution 
to simple algebraic manipulations consisting of arbitrary but fixed vectors. Because the Radon transform reduces the 
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scalar Helmholtz equation to that of simple harmonic oscillations. The poloidal, toroidal representations also suggest 
the use of the Moses eigenbasis. This simplification also enables us to write integral representations as we have seen on 
a simple example. One can try imposing different boundary conditions and symmetries motivated by considerations 
in the physical space. The CK solution is also an eigenvector of the RBS operator. 

The Euclidean topologically massive Abelian gauge fields and also the force-free magnetic fields are examples of 
Trkalian fields. Hence the Radon and the spherical curl transform analysis are also valid for these. In the topologically 
massive case, the gauge potential is also Trkalian. The Radon transform of an anti-self-dual potential (or field) is 
related by antipodal map on the sphere to the transform of the self-dual potential obtained by inverting space 
coordinates. 

The L and CK fields are also solutions of this theory. This has been overlooked previously. Furthermore, the L 
solution provides an example for quantization of the topological mass in this theory. This suggests a natural geometric 
interpretation for the Lundquist field in this context. However the discussion of physical implications of this is beyond 
our goals here. 
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APPENDIX A: THE RADON TRANSFORMATION 
1. The Radon transformation 

The Radon transform of vector fields is defined^^ as 

F'^ip, k) ^ n[F{x)\{j), i^)= j F{x)5(j) - K ■ xjd^x. (Al) 
This is a componentwise generalization of Radon transform 

/^(p, k) - n[f{x)]{p, k) I f{x)S{p - K ■ x)d'x, (A2) 
of scalar functions in cartesian coordinates. This satisfies the identities 



f'^{-p,-n) = f'^{p,K), = (A3) 

The Radon transform (|A1[) of vectors also satisfies similar identities. 

The Radon transform is simply defined^^s^ for (vector) functions in Schwartz class iS[M^] of rapidly decreasing 
functions on R'^. Hence the function vanishes at infinity as we assume repeatedly, for example in Sections III Bl and 
mil Then the Radon transform is a function in iS[S^ X M]. A natural domain of definition for the Radon transform 
is the set of hyperplanes in R'^, i. e. the projective space. The space X M is a two- fold covering of thisi^i^!^^ 
Note that the definition (|Aip actually provides an extension22i2^ of the Radon transform as a function homogeneous 
of degree --1 on (R^ — |0|) X R. The adjoint transform, see Definition [31 works in the opposite direction. We refer 
the reader to Refs. [TUB-I^ or |25| for a derivation of ([M)) . The inverse Radon transform is given as 



F{x)^--^Vl[ F'^{K-x,K)dn, (A4) 
Svr-^ Js2 

where is the Laplace operator. 

We can prove the Fourier slice theorem (|^5|) using (jAl[) . The relation (|46p for Fourier transform of the Riesz 
potential is valid for functions in the Schwartz classi^ 

The Radon probe transformation^^ is defined as 



DI[F{x)]{p,k) = V{p,k) ■1Z[F{x)]{p,k), 



(A5) 
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where V{p, n) is the probe. 

Further note that, if we had defined the Radon transformation on the projective space of hyperplanes^ in M.^, then 
in Section |V] we would take the Radon transforms of self-dual and anti- self- dual fields on distinct spaces. 

2. Inversion of the refined transformation 

The Radon transform (|A1[) is a right inverse: TZ'j^^7Z.[F]{x) — F{x) of the refinedi^ (inverse) transformation 



7^^i[F^(p, k)]{x) = -^V^ F'^iK ■ X, K 



)dn (A6) 



IH 

for any canonical hemisphere H. However as a left inverse, we find 



= -^V2 / / F'^{p,K)Sip- K-x)dndp, 



nnH^[F{p',K')]{p,K) ^ - J J F{p',K')w{p,K;p',K')dn'dp' (a?) 



F{p,k), keH, 

f{-p,-k), KeH', 



where 



W{p, p\ k') = ^ j 5{p-K- x)Vl5{p' - k' ■ x)d''x (A8) 

= -5{p' - p)5{k' - k) - 5{p' +p)5{k' + k). 

We refer the reader to Ref. [ill for evaluation of this integral. Thus only in the case F{p, k) is the Radon transform 
of a function can we define both right and left inverse of this refined transformatioujiS. 

APPENDIX B: THE SPHERICAL CURL AND RADON TRANSFORMS 

The inverse Radon transform in differential forms is given as 

^^F{x)^^A^[ *''F'^(K-x,K)dn, (Bl) 

where A^; = (d + (5)^ = 5d + d5 is the Laplace-Beltrami operator and 5ijJr = (—1)'^ * d * for an r-form aj^. 

We can prove equivalence of the spherical curl transform and the Radon transform as follows. iS, If we substitute 
([T7| in (HH), we find 



1 

We can write this as 



*«F(a;) = ^-^1/2 e-'^''^P-'^-'''>uj''{K)F^{p,K)dn. (B2) 



(2^) 

We find 



(B3) 
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following a similar reasoning given in Ref. [13 with our notation. Here we use the completeness relation ([5]) for the 
coframe, the expression ([T51) for the Radon transform of the field which yields < lu^, >= and the fact that 

the integrals of the two remaining terms over the sphere are equal. Then we are led to 

^^F{x) = J e-*^''(p-'^-^) "F'^ip, K)dn. (B4) 

One recovers the inverse Radon transform formula (jBip introducing the Laplace-Beltrami operator and restricting 
p = K ■ xJ^ 

Also substituting ([T^ tor p = k ■ x in (jBip 



* "^(^) = 2(27rV/2 5 / [^'^"'"'^''('*)^-(^^«) + e-'^'"'-^uj-{-K)sa{~Xi^K)] dn, (B5) 

we see that the inverse Radon transform (jBl[) yields the expression ((T6| of *"F(a;) in terms of the spherical curl 
transformiiS Here the second integral reduces to the first one as k — > —k. 

APPENDIX C: BIOT-SAVART INTEGRAL OPERATOR 

i) V^-BS[F]{x) = Q, 



V ■ BS[F]{x) = V -V X A[F] [x) (CI) 
V • / F{v) X ,^~^3 d^v 



= 0. 

ii) V, X BS[F]ix) = F{x), 

V^xBS[F]ix) = V^xV^xA[F]{x) (C2) 



-VxXVxX / - — ^-^d-^y, 
47r Jo \x - y\ 



4tt Jo \x - y\ 4tt Jo \x - y\ 

/ F{y)V^—^d\j-^vJ F{y)-Vy-^d'y. 

All Jo \x-y\ 47r Jo \x - y\ 



We can write the second integral as 

" F{y) 



\x - y\ 



This yields 



D \x-y\ Js\x-y 



V, X BS[F]{x) = F(a;) + -^V, / Zj^lZi^ ^3 _ 1 / Jiv) . ^ (^4) 

47r 7^, |a;-y| 47r Js\x-y\ 



47r7^ ^ ^ '\x^y\i 47r/s|a;-y|3 
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for X G D. Thus we conclude 



X BS[F]{x) = F{x), xeD, (C5) 

if 

a) F{y) is divergence-free: Vj, • F{y) — 0, 

b) F{y) vanishes at infinity, as the region D is extended to whole space [for a finite region D, F{y) is tangent 
to the surface S = dD bounding this region]. 



APPENDIX D: BIOT-SAVART INTEGRAL OF THE LUNDQUIST FIELD 

We shall evaluate 



bs[Fl]{x) = ^ / FUy) X T^^^d' 

An Jjj \x-y\-^ 



y 



(Dl) 



= h 



I A 



which can be decomposed into five terms. We choose: y = rer{(t)) + ze^ and x = Rer{9) since Fl is invariant along 
the z-axis. We shall need 



-dz = 0, 



_ (a2 + z2)3/2 

We shall also make use of Poisson integral formulas 



(„2 ^2)3/2 



dz = 2 — . 



(D2) 



1 f^'' 1 f^'' 

— P{R,r,(l}-e)d(j) = l, — Rsm(j)P{R,r. 
27r Js=o 27r J^^g 



1 

2^ 



2ir 



Rcostj) P{R,r, 



for the region r < R. Here P{R, r, i 



i<P=o 

is the Poisson kernel 



= rsin6', 
— r cos 9, 



(D3) 



P{R, r,(l)-e) 



i?2 



R^ + - 2rRcos{(l>-e)' 



(D4) 



We interchange i? and r in (ID3p for the region R < r. One can easily prove these, for example following Ref. l56|. 
For each term in (jDl[) . we shall first carry out integration over z and decompose the resulting integral into two pieces 
corresponding to regions r < R and R < r. Then we shall carry out integrations over (p using ()D3|) . We shall end up 
with integrals over r (x — vr) combining these again. The first integral 



h = - — 



zF^{r)er{(f>) 



4tt J [i?2 + r2 - 2rR cos{(j) -6) + z2]3/2 



rdrd(f)dz, 



(D5) 







immediately vanishes upon integration over z. The second term yields 



1 

An 

V 



RF,{r)eg 



[i?2 



2ri?cos( 



X 



X 



x=0 



■xdx + X 



^2J3/2 



rdrd(j)dz 
xdx 



(D6) 



x^ 



X2 



eg. 



The third integral is 
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h = — 



The fourth integral yields 



The fifth integral is 



Thus we find 



47r j [i?2 + r2 - 2ri?cos(0 -9) + z2]3/2 



rdrd(/)dz 



1 I 

-Fq— Ja{x)xdxeg - h- 



RF^{r) cos(0-6')( 



1 

47r 7 + r2 - 2r_R cos 



-If. 



22J3/2 

° Ji(a:) 



rdrdcjjdz 



x=X ^'^ 



dx 



47r J [i?2 + r2 - 2rRcos{(j) -9) + z2]3/2 
-Fq / Ji{x)dxez ~ I4. 

V Jx=X 



rdrd<j)dz 



(D7) 



(D8) 



(D9) 



— / Jo{x)xdxe0 + / Ji{x)dxe 

^ Jx=0 Jx=X 



(DIO) 



adding up these terms. Here, we use 



-j^ pX POO 

— I Jo{x)xdx — Ji{X), / Ji{x)dx — Jo 

^ Jx=0 Jx=X 



(X). 
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